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Abstract 


The  need  for  geophysical  and  astrophysical  information  is  critical  for  the 
design  of  aircraft,  missiles,  and  satellites.  The  HANDBOOK  OF  GEOPHYSICS 
AND  SPACE  ENVIRONMENTS  is  an  attempt  by  the  U.  S.  Air  Force  to  organize 
some  of  these  data  into  compact  form. 

The  effects  of  the  lower  atmosphere  on  propagation  of  electromagnetic  waves 
in  the  optical,  microwave,  and  radio  regions  are  discussed  in  this  Chapter;  the 
emphasis  is  on  information  useful  in  meteorological  investigations. 
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REFRACTION,  ATTENUATION,  AND 
BACKSCATTERING  OF  ELECTROMAGNETIC 
WAVES  IN  THE  TROPOSPHERE: 
A  Revision  of  Chapter  9,  Handbook 
of  Geophysics  and  Space  Environments 


9.1  REFRACTION  IN  THE  LOWER  TROPOSPHERE  by  V.J.  Falcone,  Jr. 

The  speed  of  propagation  of  an  electromagnetic  wave  in  free  space  is  a 

g 

constant,  c,  which  is  equal  to  3x10  m/sec.  In  a  material  medium,  such  as  the 
atmosphere,  the  speed  of  propagation  varies.  Even  small  variations  in  speed 
produce  marked  changes  in  the  direction  of  propagation,  that  is  refraction. 

In  the  atmosphere,  the  speed  of  propagation  varies  with  changes  in  composi¬ 
tion,  temperature,  and  pressure.  At  radio  wavelengths  speed  does  not  vary 
significantly  with  the  wavelength,  but  in  the  optical  region  the  speed  depends 
strongly  on  the  wavelength.  In  the  lower  15  km  of  the  atmosphere  water  vapor 
is  the  major  contributor  to  variation  in  atmospheric  composition,  and  at  radio 
wavelengths  the  speed  of  propagation  is  strongly  affected  by  water  vapor. 
Temperature  and  pressure  variations  are  principally  functions  of  altitude,  al¬ 
though  for  propagation  at  small  elevation  angles  significant  variations  may  occur 
along  horizontal  distances.  From  the  standpoint  of  effect  on  the  speed  of  propa¬ 
gation,  temperature  variations  at  any  given  altitude  are  more  significant  than 
pressure  variations. 


(Received  for  publication  14  November  1969) 


2 


In  its  most  general  form  the  refractive  index  is  a  complex  function.  The 
real  term  of  the  complex  function  is  called  the  phase  refractive  index,  n; 


c 

n  =  — 
v 


(9-1) 


where  c  is  the  speed  of  light  in  a  vacuum  and  v,  the  phase  velocity,  is  the  speed 
of  propagation  in  a  particular  medium.  In  the  troposphere  where  n  is  nearly 
equal  to  one,  it  is  convenient  to  define  the  quantity 


N  =  (n  -  1)  xlO  . 


(9-2) 


N  is  called  the  refractive  modulus;  units  of  (n  -  1)  x  10  are  called  N -units. 
9.1.1  Optical  Wavelengths 

An  approximate  relation  between  the  optical  refractive  modulus  and  atmos¬ 
pheric  pressure  and  temperature  is 


Nm  *77.6f  . 


(9-3) 


where  N  is  the  refractive  modulus  for  wavelengths  >  >  20jjl,  P  is  atmospheric 
pressure  in  millibars,  and  T  is  atmospheric  temperature  in  degrees  Kelvin. 

The  dispersion  formula  of  Edlen  (1953),  which  has  been  adopted  by  the  Joint 
Commission  for  Spectroscopy,  is 


N  ,  64.  328  +  2949^-P  +  -255a" 

s  146  -  1/  V  •  41  -  l/\ 


(9-4) 


where  Ng  is  the  refractive  modulus  at  a  wavelength  X.  for  a  temperature  of  288  K 
and  a  pressure  of  1013.  25  mb,  and  X  is  the  wavelength  in  microns.  A  somewhat 
less  precise  but  more  convenient  dispersion  formula  is 


N  =  N„ 


7.  52  x  10 
K2 


-3 


(9-5) 


Equations  9-3  and  9-5  may  be  combined  to  give  the  refractive  modulus  as  a 
functiori  of  pressure,  temperature,  and  wavelength; 


„  77.  6  P  0.  584  P 

JN  =  q,  +  O 

1  T\ 


(9-6) 
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Refractive  moduli  calculated  by  using  Eq.  9-6  will  be  in  error  no  more  than 
one  N-unit  over  the  temperature  range  -30°C  to  +30°C  for  wavelengths  from  0.  2p 
to  20  p.  Thus  Eq.  9-6  covers  the  spectrum  from  the  near  ultraviolet  through  the 
near  infrared.  Figure  9-1  is  a  nomograph  based  on  Eq.  9-6  that  gives  values  of 
N  accurate  to  about  5  N -units. 
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Figure  9-1,  Nomograph  for  Computing  Refractive  Modulus  at  Optical  Wavelengths 


9.1,2  Hadio  Wavelengths 

At  radio  wavelengths  the  relationship  of  refractive  modulus  to  pressure, 
temperature,  and  water-vapor  pressure  is  given  by: 


77.  6 
=  T  + 


373000 P 
_ wv 


(9-7) 


where  Pt  is  the  partial  pressure  of  water  vapor  in  millibars,  P  is  pressure  in 
millibars,  and  T  is  temperature  in  degrees  Kelvin.  Figures  9-2  and  9-3  are 
nomographs  based  on  Eq.  9-7;  they  give  values  accurate  to  within  5  N-units. 
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Figure  9-2.  Nomograph  for  Computing  Contribution  to  Microwave 
Refractive  Modulus  Due  to  Density  of  Atmospheric  Gas 


Figure  9-3.  Nomograph  for  Computing  Contribution  to  Microwave 
Refractive  Modulus  Due  to  Water  Vapor 
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Equation  9-7  Is  accurate  to:  0. 1  N -units  from  the  longest  radio  wavelengths  in 
use  down  to  about  6  mm  (50  GHz);  5  N-units  from  6  mm  to  4  mm  (50  to  75  GHz); 
and  1  N-unit  from  4  mm  to  2.  6  mm  (75  to  115  GHz). 

Absorption  by  atmospheric  constituents  begins  to  rise  to  significant  proportion 
with  decreasing  wavelength  beyond  1.  5  cm.  Water  vapor  content  is  by  far  the 
leading  factor  in  causing  changes  in  N,  followed  in  order  of  importance  by  tempera¬ 
ture.  and  pressure.  For  example,  for  a  temperature  of  15°C,  pressure  of  1013  mb 
near  ground  level,  and  a  relative  humidity  of  60  percent  (Pwv  =  10  mb),  the  fluctua¬ 
tion,  AN,  is 


AN  =  4.  5A  Pwv  -  1.26A  T  +  0.27A  p  (9-8) 

As  Eq.  9-8  shows,  a  fluctuation  in  water  vapor  pressure  has  16  times  the  effect 
on  the  refractive  index  as  the  same  amount  of  fluctuation  in  total  pressure,  and 
3.  5  times  the  effect  as  the  same  fluctuation  in  temperature.  Equation  9-7  may 
also  be  used  in  the  windows  of  relative  transparency  for  submillimeter  waves 
(X.  >100(1,  f  <3xiO^GHz)  with  an  error  of  10  to  20  N-units. 

9. 1.3  Standard  Profiles  of  Refractive  Modulus 

The  vertical  distribution  of  the  refractive  modulus  can  be  calculated  from 
Eq.  9-7  using  vertical  distributions  of  vapor  pressure  and  temperature  as  a 
function  of  pressure. 

Under  normal  conditions,  N  tends  to  decrease  exponentially  with  height.  An 
exponential  decrease  is  usually  an  accurate  description  for  heights  greater  than 
10,000  feet  (3  km);  below  10,000  feet,  N  may  depart  considerably  from  exponential 

I 

behavior.  The  median  value  for  the  gradient  dN  is  typically  -  0.  012  N-units  per 
foot  for  the  first  few  thousand  feet  above  ground  level. 

For  many  purposes  it  is  desirable  to  have  standard  refractive -modulus  profiles 
for  the  atmosphere.  By  using  the  equations  of  the  model  atmosphere  (Chapter  2), 
an  exact  analytical  expression  for  the  standard  optical  refractive  modulus  can  be 
derived.  A  simplified  approximation  to  this  is 

Nqo  =  273  exp(-Z/  32.  2),  (Z<25);  (9-9) 

Z  is  the  altitude  in  thousands  of  feet. 

Equation  9-9  can  be  differentiated  to  obtain  the  standard  gradient  of  optical 
refractive  modulus; 

dN 

“T7~=  -8.  48  exp(-Z/  32,  2),  (Z<  25)  .  (9-10) 
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Equations  9-9  and  9-10  may  be  corrected  for  dispersion  through  use  of  Eq.  9-5. 

For  the  radio  wavelengths  it  is  necessary  to  assume  a  distribution  of  water 
vapor  in  order  to  obtain  an  expression  for  the  refractive  modulus.  Assuming 
Pwv  m  10.  2  (1  -  0. 0195  Z)^,  for  Z  <  25,  a  simplified  approximation  is 

N  =  316  exp(-Z/26.  5),  (Z  <  25)  .  (9-11) 

The  standard  gradient  of  radio -wave  refractive  modulus  is  then: 

-jj§=  -11.9  exp(-Z/26.  5),  (Z<25).  (9-12) 

Figures  9-4  and  9-5  are  graphs  of  standard  profiles  calculated  from  Eq.  9-10 
through  9-12. 

9. 1.4  V'ariations  of  Refractive  Moduli 

Actual  profiles  may  differ  markedly  from  the  standard  profiles.  Figures  9-6 
through  9-9  show  some  profiles  of  refractive  modulus  at  microwave  frequencies 
calculated  from  radiosonde  measurements.  These  are  considered  typical  for  the 
air  masses  indicated.  Average  deviations  from  a  model  atmosphere  refractive 
index  have  been  studied  extensively:  for  example,  see  Bean  and  Dutton  (1968). 

Cumulus  clouds  are  evidence  of  the  existence  of  a  very  inhomogeneous  atmos¬ 
phere.  Figure  9-10  shows  some  measurements  of  refractive  modulus  and 
associated  parameters  within  a  fair-weather  cumulus  cloud.  The  time  response 
of  the  instruments  from  which  refractive  modulus,  temperature,  and  water-vapor 
pressure  were  obtained  was  such  that  changes  occurring  in  distances  as  small  as 
5  feet  (1.  5m)  could  be  measured.  However,  the  instrument  for  measuring  liquid 
water  content  was  much  slower  responding.  Figure  9-11  shows  a  composite  cloud 
which  summarizes  data  from  30  cloud  passes. 

Figure  9-12  shows  the  average  An  cloud  to  clear  air  to  be  expected  in  various 
parts  of  the  United  States  at  the  midseason  months.  The  chances  of  having  cumulus 
clouds  at  1500  h  local  time  for  these  months  is  also  shown.  Additional  climatologi¬ 
cal  data  on  AN  and  cumulus  clouds  is  given  by  Cunningham  (1962). 

The  deviations  in  refractive  modulus  are  principally  in  the  vertical  direction. 
Regions  of  more  or  less  constant  gradient  of  the  refractive  modulus  are  called 
stratified  layers.  The  horizontal  extent  of  these  layers  may  vary  from  a  few  miles 
to  hundreds  of  miles  depending  on  the  meteorological  processes  by  which  they  are 
produced.  When  N  decreases  with  height  inside  a  specific  layer  much  faster  than 
it  does  above  or  below  the  layer,  the  layer  is  said  to  be  super -refracting  for 


propagation.  One  cause  of  this  is  temperature  inversion  (Section  3.  1.5).  Layers 
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Figure  9-5.  Variation  of 
Standard  Refractive  Modulus 
with  Altitude 


100  140  ISO  220  260  300  340  380 

MODULUS  (N  units) 


Figure  9-4.  Variation  of  Standard 
Gradient  of  Refractive  Modulus  with 
Altitude 
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Figure  9-6.  Microwave  Refractive 
Modulus  Profile  in  Continental  Tropical 
Air  Mass 
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Figure  9-7.  Microwave  Refractive 
Modulus  Profile  in  Continental 
Polar  Air  Mass 


Figure  9-8.  Microwave 
Refractive  Modulus  Profile 
in  Continental  Polar  Air  Mass 
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Figure  9-9.  Microwave  Refractive 
Modulus  Profile  in  Maritime  Polar 
Air  Mass 
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Figure  9-10.  Aircraft  Measurements  Through  a  Cumulus  Cloud.  The  heavy  lines 
show  the  time  during  which  an  observer  in  the  aircraft  indicated  that  the  plane  was 
within  the  visible  cloud.  The  calculated  virtual  temperature,  Tv,  was  corrected 
for  liquid  water  content;  T  is  the  measured  temperature.  Relative  humidity  in 
percent  is  shown  on  the  ambient  water-vapor  pressure  curve;  the  curve  labelled 
es  is  the  calculated  saturation  water-vapor  pressure  at  the  temperature  encountered. 
The  bottom  curve  shows  the  liquid  water  content  (LWC) 
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Figure  9-12  (Cont. )  Percent  Frequency  of  Cumulus  and  Cumulonibus 
Cloud  (Dashed  Contours)  and  A  N,  the  Average  Change  in  Refractive 
Modulus  (Solid  Contours)  for  January,  April,  July  and  October 
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of  negative  height -gradient  of  N  In  association  with  regions  in  which  the  tempera¬ 
ture  gradient  is  positive  (or  less  negative)  than  the  gradients  of  the  layers  just 
above  and  below  are  known  as  subsidence  inversion  layers.  These  layers  generally 
have  a  large  horizontal  extent.  Layers  in  close  proximity  to  the  earth's  surface 
are  strongly  influenced  by  the  local  conditions  of  the  earth's  surface  and  for  this 
reason  are  less  extensive  than  the  layers  described  above. 

In  the  inversion  layer,  the  temperature  may  change  by  a  few  degrees  in  in¬ 
tervals  of  from  fifty  to  a  few  hundred  feet  in  altitude.  This  temperature  change  ' 
accounts  for  only  a  small  change  of  a  few  N -units.  However,  an  inversion  usually 
indicates  the  presence  of  a  humid  air  mass  under  a  dry  one.  The  transition  from 
humid  to  dry  air  causes  a  marked  change  of  N  in  the  super-refracting  layer, 
typically  of  20  to  50  N-units;  changes  as  pronounced  as  80  N-units  have  been 
measured.  Super -refracting  layers  may  be  clear -weather  phenomena,  or  can  be 
accompanied  by  haze  (aerosols)  in  the  lower  air  mass.  Invariably  they  signify 
stable  weather  situations,  such.as  occur  when  a  high  pressure  center  stagnates  in 
an  area. 

The  horizontal  and  temporal  extent  of  super-refractive  layers  varies  widely. 

In  New  England  it  may  be  only  a  few  tens  of  miles.  In  midwestern  states,  the  layers 
extend  farther  and  may  last  from  a  half  hour  up  to  a  week.  "In  the  trade-wind  zones 
of  the  world,  the  climatic  regime  (manifested  by  steady  wind  directions  and  speeds 
throughout  most  of  the  year)  sustains  super -refracting  layers  which  extend  a  few 
thousand  miles  both  east  to  west  and  north  to  south. 

9.1.5  Turbulence 

It  has  long  been  recognized  that  the  effects  produced  by  the  atmosphere  upon 
electromagnetic  waves  propagating  through  it  are  a  measure  of  the  nature  of  the 
atmosphere.  If  the  atmosphere  is  considered  as  a  medium  with  electromagnetic 
properties  that  are  random  functions  of  space  and  time,  then  atmospheric  properties 
can  be  investigated  by  measurements  of  wave  propagation;  these  observations  are 
called  remote  probing.  The  atmospheric  (or  meteorological)  parameters  are  ' 
inferred  from  their  influence  on  Eqs.  9-6  and  9-7;  it  is  assumed  that  scattering  is 
due  to  random  fluctuations  in  the  dielectric  constant  of  the  atmosphere.  When 
fluctuations  in  the  refractivity  are  of  interest  they  are  studied  via  the  correlation 
function  and  its  Fourier  transform,  the  spectral  density.  This  approach  is  described 
by  Tatartski  (1961  and  1967),  Staras  and  Wheelon  (1959),  Hufnagle  and  Stanley  (1964), 
and  Strohnbehn  (1968). 

A  locally -homogeneous  isotropic  turbulent  model  of  the  atmosphere  is  assumed, 
that  is  a  model  of  well  mixed  random  fluctuations.  This  model  is  restrictive  and 
requires  justification  for  each  new  application.  One  assumes  the  spectral  density 
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$  n  (  k  )  (the  three  dimensional  spectrum  of  refraciivity)  is  given  by 

$n(K  )=  0.033  C2*  'H/3exp^-K2/Km2  j  [cm3]  (9-13) 

where  K0  <  K  <Km  •  The  spatial  wave  number,  k  (cm-1)  is  2ir/f  ,  f  is  the  size  of 

the  turbulent  eddy  or  blob,  k  is  Zir /  L  ,  L  is  the  outer  scale  of  turbulence  (typically 

3  4  °  °  0  -1 

the  order  of  10  to  10  cm,  depending  on  layer  height),  Km  is  5.  92 /f  (cm  ),  f 

is  the  inner  scale  or  turbulence,  and  CR  is  specified  by  Eq.  9-14.  Figure  9-13 
shows  the  typical  spectrum  of  irregularities,  #„(•<),  and  the  ranges  of  energy 
input,  redistribution,  and  dissipation.  Physically,  the  energy  is  put  into  the  tur¬ 
bulence  from  the  largest  scale  sizes  (smallest  value  of  k  )  by  wind  shear  dnd  .con  - 
vective  heating;  the  energy-producing  eddies  are  assumed  to  have  a  spacial  wave 
number  less  than  kq.  The  region  between  *  and  K  m  is  the  redistribution  (inertial) 
range,  where  energy  is  transferred  from  large  eddies  (small  k)  to  smaller  eddies 
(larger  k  )  until  viscous  effects  becomes  important  at  k  =  5.  92/ 1  and  the  energy 

is  dissipated.  Near  the  ground,  i  is  the  order  of  0.  1  to  1  cm. 

2  o  2 

Once  the  value  of  Cn  is  determined,  the  spectral  density  is  known.  Cr  may  be 

found  directly  from  the  dimensionless  structure  function  Dn(r), 


Dn(#)  =  [n  (?+?j)  -  N  (rj)]2  =  C2  r2/3  (9-14) 

where  N  (r)  is  the  normalized  fluctuating  part  of  the  index  of  refraction,  the  bar 
indicates  the  average  of  the  squared  quantity,  and  r  (r  >>  f  )  is  the  size  of  the  in- 
homogenities  determined  from  the  difference  in  the  values  of  N  at  two  points  ?  and 


At  optical  wavelengths,  Cn  is  determined  from  temperature  measurements 
alone; 

Cn  =  10'6  (p/p0)CT  [cm-1/3]  (9-15) 

where  p/pQ  is  the  ratio  of  the  average  atmospheric  density  at  a  given  altitude  to 
the  density  at  sea  level.  The  structure  constant,  Cj,,  is 

CT  =  2.4  e1/3  (y/p)  [°Ccm'1/3]  (9-16) 

2  -3 

where  e(cm  )  is  the  rate  of  energy  per  unit  mass  dissipated  by  viscous  friction, 
y(°C  cm*?)  is  the  average  vertical  gradient  of  the  potential  temperature,  and 
P  (s  *)  is  the  average  shear  rate  of  the  wind  (Haufnagle  and  Stanley,  1964). 


15 


Figure  9-13.  Three-dimensional  Spectrum  of  Refractive  Index 
Fluctuations.  Note  that  units  for  the  ordinate  scale  are  m11/  3. 
(After  Strohnbehn,  1968.) 


Figure  9-14.  Dissipation  Rate  e  vs  Altitude.  The  curve  corresponds 
to  an  estimated  average  e.  (After  Hufnagel  and  Stanley,  1964.) 


\ 
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Figure  9-14  show*  the  average  dissipation  rate  estimated  from  observed  values 
as  a  function  of  altitude.  Figures  9-15  and  9-16  are  plots  of  observed  geometric- 
mean  values  of  y  and  p  respectively  as  functions  of  altitude.  The  values  of  y  are 
taken  from  the  1966  Supplementary  Atmosphere;  the  values  of  f3  are  computed  from 
reported  wind  profiles.  Figure  9-17  shows  CR  as  a  function  of  altitude;  the  values 
of  p/pQ  are  taken  from  the  1966  Supplementary  Atmosphere.  At  radio  wavelengths, 
Eq.  9-15  cannot  be  used  because  water  vapor  must  be  considered  as  well  as 
temperature  (Crane,  1968). 

Cn  for  a  model  atmosphere  can  be  obtained  directly  from  Figure  9-17.  When 
values  of  p/  pQ  are  known  from  radiosonde  observations,  CR  for  the  given  p/  pQ  can 
be  calculated  by  using  values  obtained  from  Figures  9-14  through  9-16  in  Eqs.  9-15 
and  9-16.  The  most  accurate  determination  of  Cr,  however,  is  obtained  from  back- 
scatter  measurements. 

Tropospheric  scatter  experiments  are  used  to  probe  the  atmosphere,  see 
Staras  und  Wheelon  (1959)  and  Hardy  and  Katz  (1969).  At  radio  frequencies,  the 
mechanism  responsible  for  backscatter  and  forward  scatter  .beyond  the  horizon  is 
the  refractive  index  variation  due  to  fluctuations  in  properties  of  the  atmosphere. 
The  ratio  of  the  received-to-transmitted  power  (P  J  Pt)  depends  upon  the  integral 
of  the  scattering  cross  section  per  unit  volume  over  the  common  volume  defined 
by  the  antenna  pattern  or  patterns  for  backscatter  and  forward  scatter  respectively. 

Pr/Pt  .  (X2/ 4w2)  f  (G^^rJr^)  pd3r  (9-17) 

Volume 

where  X  is  the  wavelength  employed,  and  Gr  are  the  gains  of  the  transrhitting 
and  receiving  antennas  respectively,  Rj  and  Rg  are  the  distances  from  the  scatter¬ 
ing  volume  (d3r)  to  the  respective  antennas,  and  fi  (a  reciprocal  length)  is  the 
scattering  cross  section  per  unit  volume.  Figure  9-18  illustrates  the  path  geometry 
for  Eq.  9-16. 

The  scattering  cross  section  per  unit  volume  is  directly  related  to  the  spectral 
density: 

p  =  (ir2/X4)*n(K)  (9-18) 


where  |x|  is  (4ir/X)  sin  (8/2). 


Scattering  is  not  the  only  effect  of  atmospheric  turbulence  on  the  propagation 
of  electromagnetic  waves.  As  the  waves  propagate  through  the  atmosphere, 
fluctuations  in  amplitude  and  phase  occur.  The  amplitude  and  phase  fluctuations 


.4 


v 

N 

U 


K 

►- 


1.4(10 


I.2mI0-4 


Sf  I.O*IO-« 
o 


8*I0-9  [■ 
6*10*® 
4(10“  9 


*- 
< 
K 
U J 
Q. 

a 

Ui 


-J 

< 


O 

0. 


2x10 


-5 


■i  i  i  1 1 


1  *  i  i  i  i  1 1 


Ikm 


10km 

ALTITUDE 


100km 


Figure  9-15.  Average  Potential  Temperature  Gradient  (\)  vs 
Altitude.  Computed  from  the  U.  S.  Standard  Atmosphere.  (After 
Hufnagel  and  Stanley,  1964.) 


Figure  9-16,  Average  Wind  Shear  (p)  vs  Altitude.  Computed  for 
0.  5 -km  vertical  separations.  (From  Hufnagel  and  Stanley, 

1964.) 


ALTITUDE 


Figure  9-17.  Index  of  Refraction  Structure  Constant  (Cn)  vs 
Altitude.  (From  Hufnagel  and  Stanley,  1964.) 


«V 


Figure  9-18.  Geometry  for  Turbulent  Scattering. 
Rl  and  R2  are  the  distances  from  the  scattering 
volume  d*r  to  the  transmitter  and  receiver, 
respectively,  and-©-  is  the  scattering  angle 
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may  be  described  using  the  equations  of  geometrical  optics  or  the  smooth  pertur¬ 
bation  solution  of  the  full  wave  equation,  see  Tatarskl  (19611,  a  summary  of  which 
follows. 

In  order  to  apply  the  geometrical  optics  approximation,  the  conditions  that 
have  to  be  satisfied  are:  — 

\  «fQ  and  (XL)1/2  «fD  , 

where  X  is  the  wavelength  of  the  propagating  wave,  tQ  is  the  inner  scale  of  tur¬ 
bulence,  aipd  L  is  the  path  length  over  which  the  wave  propagates.  The  phase 
fluctuation* at  a  point  r  along  the  path  is 


S  (?)  =  k  f  n(s)  ds 


(9-19) 


where  S(?)  is  the  total  phase  change,  k  is  2ir/X,  n(s)  is  the  index  of  refraction  in 
the  direction  of  t,  and  ds  is  the  element  of  path  length.  The  amplitude  fluctuation, 
x ,  at  point  ?  along  the  path  is 


X=  in  [a(?)/Aq(?)] 

x  =  kfQ  VT  S1  (  5  *  y*  z)  d5 


(9-20) 


w  tere  A  ( r )  is  total  amplitude  at  point  r,  is  the  transverse  Laplacian,  and 
£  is  the  direction  of  propagation.  Sj  (  ? ,  y,  z)  is  the  phase  fluctuation  about  its 
mean  value. 


/r 

6n(s)  ds 


(9-21) 


where  6n(s)  is  the  deviation  of  the  index  of  refraction  from  its  mean  value.  If 
n(s)  depends  only  on  altitude,  then  Eq.  19  may  be  used  to  determine  the  average 
phase  change  by  substituting  the  average  index  of  refraction  in  the  integrand. 


*  fl 


(s)  >  ds;  this  is  useful  in  studying  refraction. 


The  covariance  function  of  phase  fluctuations,  C  (  5  ),  at  two  receivers  that 

s 

are  both  an  equal  distance  L  away  from  the  transmitter  and  are  separated  from 
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each  other  by  a  distance 


/  2  2. 
(y  +  z  ) 


1/2 


is 


00 

FJK  )J  (KC)KdK 
B  O 


(9-22) 


J  (  k  5)  is  the  Bessel  function  and  F  ( k)  is  the  two  dimensional  Fourier  transforms 
of  Cfl  (  ?); 

F„(k)  .  2irk2  L  *  (k)  [cm2]  (9-23) 

B  lx 

The  covariance  function  of  amplitude  fluctuations,  Cx  (  ?)»  for  the  conditions  given 
above  is 


C*<C) 


Fx(K)Jo(KS)KdK 


and 


(9-24 


Fx(k)  »  (itA3/6)k4  ®n(K) 

[cm2) 

(9-25) 

Under  restrictions  that  X  <<  f  Q  and  that  L  <<  i  X^  , 
solution  of  the  full  wave  equation  leads  to 

the  smooth  perturbation 

Fs(k  )  »  irk2  L  j 

[cm2) 

(9-26) 

and 

Fx(k)  =  irk2  L  j 

'i-8*™  )»„(«> 

[cm2] 

(9-27) 

2 

where  a  is  k  L/k.  These  restrictions  limit  the  validity  of  Eqs.  9-26  and  9-27 
to. wavelengths  in  the  millimeter  and  optical  regions  and  to  relatively  short  paths, 
although  in  certain  cases  these  restrictions  may  be  relaxed,  -  for  example  for 
X  >  f0  .  C  and  F  can  be  measured  only  over  a  finite  range  of  values  of  k  ,  there¬ 
fore  a  complete  knowledge  of  C  or  of  F  is  not  possible. 


9.2  ATTENUATION  AND  BACKSCATTERING  IN  THE  TROPOSPHERE  by  R.  Dyer 

Scattering  and  attenuation  are  in  general  complicated  functions  of  particle 
size  and  dielectric  properties.  The  square  root  of  the  dielectric  constant,  m,  is 
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V^nn-lK  (9_28) 

where  n  is  the  phase  refractive  Index  and  Kthe  absorption  index  of  the  medium, 
i  is  V“1  • 

9.2.1  Dielectric  Properties 

It  is  convenient  to  use  the  parameter  K,  defined  by 

K  *  -^4  .  (9-20) 

in  describing  the  dielectric  properties  of  the  troposphere.  When  the  wavelength 
is  long  in  comparison  with  the  size  of  the  particles  (Rayleigh  approximation),  the 
backscatter  and  absorption  are  relatively  simple  functions  of  K.  In  this  special 
case,  backscatter  is  proportional  to  and  the  attenuation  is  proportional  to 

the  imaginary  part  of  minus  K,  Im(-K). 

The  variation  of  the  dielectric  characteristics  of  water  with  temperature 
and  wavelength  have  been  tabulated  by  Gunn  and  East  (1954).  For  water,  |k|  2 
is  practically  constant;  over  a  wide  range  of  temperatures  and  wavelengths  in 
the  centimeter  range  |k|  2  =  0.  93.  Similarly,  |k|  2  ■  0.  176  for  ice  of  normal 
density  (0.  917  g/cm\ 

9. 2.2  Backscatterinp  and  Attenuation  Cross  Sections 

The  echo  power  returned  by  a  scattering  particle  is  proportional  to  its 
backscattering  cross  section,  o.  The  power  removed  by  an  attenuating  particle 
is  proportional  to  the  total  absorption  cross  section.  Q^.  The  size  parameter 
(electrical  size)  is  xD/k;  D  is  the  particle  diameter  and  X.  the  wavelength  of  the 
incident  radiation.  When  the  diameter  of  the  scattering  or  attenuating  particle  is 
small  with  respect  to  k,  the  backscattering  and  total  absorption  cross  sections 
may  be  expressed  with  sufficient  accuracy  by  the  Rayleigh  approximation  (molecular 
scattering,  see  Section  7.  4). 

For  spherical  particles,  if  xD/k  <  0.  2, 


IkI 


V 


[cm2] 


(9-30) 


For  particles  with  ffD/k  >0.2,  a  should  be  computed  from  the  equations  of  the  Mie 

theory  of  scattering  (Battan,  1959).  Figures  9-19  and  9-20  show  normalized  back- 

2 

scattering  cross  section  (4a/ xD  )  for  ice  and  for  water  versus  the  size  parameter. 
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computed  from  the  exact  Mie  equations.  The  normalized  curve  for  ice  is  invariant 
with  wavelength  in  the  microwave  region;  the  normalized  curve  for  water  is  for 
water  at  0°C  at  3.  2 -cm  wavelength.  As  Figures  9-19  and  9-20  show,  ice  spheres 
equal  to  or  larger  than  the  wavelength  may  scatter  more  than  an  order  of  magnitude 
greater  than  water  spheres  of  the  same  size;  this  is  confirmed  by  experimental 
measurements  (Atlas  et  al. ,  1960). 

Measurements  at  5  cm  wavelength  indicate  that  the  backscattering  of  so-called 
"spongy"  hail  (a  mixture  of  ice  and  water)  is  3  to  4  dB  above  that  of  the  equivalent 
all-water  spheres  and  at  least  10  dB  above  that  of  the  equivalent  solid  ice  spheres 
(Atlas  et  al.  ,  1964).  Because  of  the  variabilities  of  sizes,  shapes,  and  liquid 
water  content  of  hail,  no  general  rules  concerning  backscattering  and  attenuation 
cross  sections  for  hail  can  be  made.  As  a  first  approximation,  however,  the  ice 
curve  of  Figures  9-  19  and  9-20  may  be  used. 

For  spherical  particles,  if  ffD/X.  <  0.  1, 


Q* 


*2d3 


lm 


(-K)+-|o  [cm2] 


(9-31) 


When  *D/\  >  0.  1,  Q.  must  also  be  computed  from  the  exact  Mie  equations.  Figure 
9-21  presents  curves  of  the  normalized  total  attenuation  cross  section  (4Qt/»rD  ) 
for  water  at  0°C  at  two  wavelengths.  Figure  9-22  shows  the  normalized  total 
attenuation  cross  section  for  ice. 


9.2.3  Reflectivity 

The  average  echo  power  returned  by  a  group  of  randomly  distributed  scatter¬ 
ing  particles  is  proportional  to  their  reflectivity,  r).  Reflectivity  is  defined  as  the 
summation  of  the  backscatter  cross  sections  of  the  particles  over  a  unit  volume; 
rj  3  So  .  When  the  backscattering  particles  are  spheres  and  are  small  enough' 
with  respect  to  wavelength  so  that  the  Rayleigh  approximation  can  be  used  (that 
is  jtD/\<  0.2),  the  reflectivity  is  proportional  to  the  radar  reflectivity  factor,  Z, 
which  is  the  summation  over  a  unit  volume  of  the  sixth  power  of  the  particle 

e 

diameters;  Z  s  SD  .  Summation  over  a  unit  volume  of  Eq.  (9-30)  gives. 


Sa  ■  q 


x  10 


-12 


[cm'1) 


(9-32) 


6  _  3 

for  Z  in  conventional  units  of  mm  m  and  *  in  centimeters. 

When  the  particles  are  larger  than  Rayleigh  size  or  composed  of  ice  or  water- 
ice  mixtures,  it  is  common  practice  to  measure  the  radar  reflectivity  and  express 
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Figure  9-19.  Calculated  Values  of  the  Normalized  (4a/irD^)  Backscatter  Cross 
Section  for  Water  at  3.  21  cm  and  0°C  and  for  Ice  at  Wavelengths  from  1  to  10  cm 


Figure  9-20.  Calculated  Values  of  the  Normalized  (4o/jtD^)  Backscatter  Cross 
Section  for  Water  at  3.  21  cm  and  0°C  and  for  Ice  at  Wavelengths  from  1  to  10  cm 
(Detail  of  Figure  9-19. ) 


*D/X 

2 

Figure  9-21.  Normalized  (4Qt/?rD  Attenuation  Cross  Sections  of 
Water  Spheres  at  0°C  and  Wavelengths  of  3.  2  and  5.  5  cm.  (After 
B.  M.  Herman,  S.  R.  Browning  and  L.  J.  Battan,  Tech.  Rept. 

No.  9,  Institute  of  Atmospheric  Physics,  University  of  Arizona, 
1961.) 
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n 

Figure  9-22.  Normalized  (4Qt/?rD  Attenuation  Cross  Sections  of 
Ice  Spheres  at  Wavelengths  Between  1  and  10  cm.  (After  B.  M. 
Herman  and  L.  J.  Battan,  Proc.  Ninth  Weather  Radar  Conference, 
Amer.  Meteorological  Society,  p.  259,  1961.) 
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it  in  terms  of  an  equivalent  reflectivity  factor,  Zg.  By  substituting  J K | 2  ■  0.  93 
(water  in  the  centimeter  wavelength  range)  in  Eq.  (9-32). 

Zg  »  3.5  x  10^  (mm®m  .  (9-33) 

Thus  Zg  is  simply  2  D  of  water  drops  of  Rayleigh  size  that  would  have  the  same 
reflectivity  as  that  measured. 

Because  Z  is  a  meteorological  parameter  that  depends  only  on  the  particle 
size  distribution  and  concentration,  it  is  useful  to  express  Z  as  a  function  of  the 
precipitation  rate,  R.  Table  5-5  (Section  5.2  Chapter  5  Revised,  Air  Force 
Surveys  in  Geophysics  No.  212  ,  AFCRL-69-0487  ,  October  1969)  lists  Z 

as  a  function  of  R  for  various  types  of  rain  at  various  locations  in  the  Northern 
Hemisphere. 

The  scattering  properties  of  snow  are  complex.  Snow  may  occur  in  a  great 
variety  of  forms  of  either  single  ice  crystals  or  aggregates  of  such  crystals. 
Measurements  suggest  the  following  relationships  are  reasonable: 

Z  =  500  R*'  6,  for  single  crystals 

and 

Z  =  2000  R^'  for  aggregates 

R  is  the  snowfall  rate  in  millimeters  of  water  per  hour. 

Carlson  (1968)  suggests  that  a  reasonable  average  for  all  types  of  snow 
measured  aloft  would  be 

Z  =  1000  R2'  67  (9-36) 


(9-34) 

(9-35) 


Clouds  composed  of  water  particles  scatter  very  poorly.  Nevertheless, 
high-power,  short-wavelength  radars  may  sometimes  detect  water  clouds  at  short 
ranges.  An  empirical  relationship  for  water  clouds  is 

Z  =  0.  048  M2,  (9-37) 

where  M,  the  water  content,  is  in  grams  per  cubic  meter. 

The  reflectivity  of  ice  clouds  is  not  known.  However,  many  so-called  ice 
clouds  are  really  snow  crystals  falling  and  evaporating  in  a  dry  layer  below  the 
cloud  base. 
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9.2.4  \ttenuntion 

Ir<  the  intensity  of  radiation  after  travelling  a  distance  r  through  an  atmos¬ 
pheric  path  of  constant  attenuation,  can  be  related  to  the  initial  intensity  I  by  a 
power  law, 

!r  8  lo  10'Ar/l°  ;  (9-38) 

A  is  the  attenuation  coefficient  in  units  of  dB  per  kilometer,  and  r  is  in  kilometers. 
A  is  the  sum  of  the  attenuation  coefficient  due  to  precipitation  (A^),  the  attenuation 
coefficient  due  to  clouds  (Ac),  and  the  attenuation  coefficient  due  to  atmospheric 
gases  (A  );  A  =  A  +  A  +  A  .  For  radar  (two-way  transmission),  the  distance  r 
is  twice  the  range  (twice  the  path  length  from  source  to  target). 

9.2.4.  1  ATTENUATION  BY  PRECIPITATION 

In  recent  years,  measurements  of  the  rainfall  attenuation  in  the  millimeter 
and  centimeter  wavelength  range  have  been  consistently  higher  than  the  attenuations 
predicted  theoretically.  Much  of  the  discrepancy,  however,  is  probably  caused  by 
an  inadequate  representation  of  the  rain  field  during  the  attenuation  experiments. 
Figure  9-23  shows  a  summary  of  the  results  of  12  such  investigations  at  10  dif¬ 
ferent  wavelengths  between  4.  3  mm  and  3.2  cm.  The  solid  lines  show  the  maximum 
and  minimum  measured  attenuation  coefficients  per  unit  rainfall  rate.  The  dashed 
line  is  the  measured  mean,  and  the  hatched  line  is  the  mean  that  would  have  been 
predicted  from  theory.  Attenuation  is  a  function  of  temperature  as  well  as  of 
wavelength  and  precipitation  rate  (or,  more  exactly,  drop  size  distribution). 
However,  the  variability  in  the  measured  attenuations  exceeds  the  temperature  cor¬ 
rection,  and  it  is  reasonable  to  assume  (Medhurst,  1965)  that  the  attenuation  co¬ 
efficient  is  a  linear  function  of  rainfall  rate.  Measurements  of  the  attenuation 
caused  by  snow  are  difficult  to  obtain  because  the  effect  is  so  small.  Figure  9-24 
shows  the  theoretical  maximum  attenuation  coefficients  assuming  a  maximum  snow¬ 
fall  rate  of  10  mm  of  water  per  hour.  Because  snowfall  rates  seldom  exceed 
3  mmwh  attenuation  due  to  snow  should  generally  be  one-third  or  less  the  value 
shown  in  Figure  9-24.  For  more  exact  computations  of  the  attenuation  to  be  ex¬ 
pected  at  various  temperatures,  wavelengths  and  snowfall  rates,  see  Gunn  and 
East  (1954). 

9.  2.  4.  2  ATTENUATION  BY  CLOUDS 

Attenuation  by  clouds  depends  on  the  liquid  water  content,  the  temperature, 
and  whether  the  cloud  is  composed  of  ice  or  water  particles.  Figure  9-25  gives 
representative  attenuations  per  unit  water  content  of  clouds  composed  entirely  of 
water  particles,  and  of  clouds  composed  of  ice  crystals.  Note  the  temperature 
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Figure  9-23.  Attenuation  (dB/km)  per  Unit  Rainfall  Rate  (mm/h) 
as  a  Function  of  Wavelength.  The  solid  lines  are  the  maximum  and 
minimum  attenuations  measured  at  each  wavelength,  the  dashed 
line  is  the  mean  of  the  measurements,  and  the  hatched  line  is  the 
mean  predicted  from  theory.  (After  Medhurst,  1965) 


and  wavelength  dependence,  and  the  sharp  decrease  in  attenuation  in  going  from 
the  liquid  to  the-solid  phase. 

9.  2.  4.  3  ATTENUATION  BY  ATMOSPHERIC  GASES 

Molecular  oxygen  and  water  vapor  are  the  only  atmospheric  gases  causing 
measurable  attenuation  of  microwaves.  Atmospheric  pressure,  temperature,  and 
water-vapor  content  are  the  governing  variables.  There  is  an  oxygen  absorption 
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Figure  9-24.  Theoretical  Attenuation  Due  to  Snowfall  of  10  mm 
of  Water  Content  per  Hour  as  a  Function  of  Wavelength  and 
Temperature.  Snowfall  rates  rarely  exceed  3  mm  of  water  content 
per  hour,  hence  usual  attenuations  should  be  approximately  one- 
third  or  less  of  these  maximum  values 


line  at  0.  5  cm  and  a  water  vapor  absorption  line  at  1.  35  cm.  Table  9-1  shows 
the  variation  of  attenuation  by  oxygen  with  wavelength  (0.  7  to  10  cm)  at  20°C  and 
one-atmosphere  pressure.  Attenuation  by  oxygen  is  pressure  and  temperature 
dependent;  it  increases  slowly  as  the  temperature  decreases,  and  increases  as 
the  square  of  the  pressure.  Table  9-2  lists  correction  factors  to  be  used  at  var¬ 
ious  temperatures.  The  attenuation  given  in  Table  9-1  must  be  multiplied  by  a 
factor  given  in  Table  9-2  in  order  to  correct  for  pressure  and  temperature  when¬ 
ever  these  variables  are  not  equal  to  20°C  and  one-atmosphere  respectively.  Both 
Tables  9-1  and  9-2  assume  a' constant  ratio  of  molecular  oxygen  to  other  atmos¬ 
pheric  gases;  this  assumption  is  valid  up  to  altitudes  of  approximately  80  km. 
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Figure  9-25.  Attenuation  (dB/km)  per  Unit  Water-content  (g/m  ) 
of  Ice  and  Water  Clouds  as  a  Function  of  Wavelength  and  Tempera¬ 
ture.  The  solid  lines  are  the  attenuation  curves  for  clouds  composed 
of  ico  particles,  and  the  dashed  curves  are  for  water  clouds 


Water  vapor  is  a  much  more  variable  constituent  of  the  atmosphere  than  oxygen. 
Table  9-3  tabulates  the  attenuation  as  a  function  of  temperature  and  wavelength 
(Gunn  and  East,  1954);  the  relation  is  only  approximate  for  operating  wave¬ 
lengths  in  the  vicinity  of  1.  35  cm  because  of  the  water  vapor  absorption  line 
at  1.  35  cm. 
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Table  9-1.  Attenuation  Due  to  Molecular  Oxygen  at  a  Temperature 
of  20°C  and  a  Pressure  of  1  -Atmosphere 


Wavelength 

(cm) 

Attenuation 
(dB  /  km) 

10 

6.  5  x  10‘3 

7.5 

7.0  x  10"3 

3.2 

7.  2  x  10~3 

1.8 

7.  5  x  10‘3 

1.5 

8.  5  x  10’3 

1.25 

1.4  x  10-2 

0.8 

7.  5  x  10'2 

0.7 

1. 9  x  10° 

Table  9-2.  Correction  Factors  for  Oxygen  Attenuation 


Temperature 

(°C) 

jb 

Correction  Factor 

20 

1.00  P2 

0 

1.19  P2 

-20 

,  1.45  P2 

-40 

1.78  P2 

* 

P  is  pressure  in  atmospheres. 
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The  effects  of  the  lower  atmosphere  on  propagation  of  electromagnetic  waves 
in  the  optical,  microwave,  and  radio  regions  are  discussed  in  this  Chapter;  the 
emphasis  is  on  information  useful  in  meteorological  investigations.  '■ 
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